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Like a silver thread, quantum entanglement [lj runs through the founda- 
tions and breakthrough applications of quantum information theory. It cannot 
arise from local operations and classical communication (LOCC) and there- 
fore represents a more intimate relationship among physical systems than we 
may encounter in the classical world. The 'nonlocal' character of entanglement 
manifests itself through a number of counterintuitive phenomena encompass- 
ing Einstein-Podolsky-Rosen paradox j2l|3j, steering [4|, Bell nonlocality [5j or 
negativity of entropy |51 [7J. Furthermore, it extends our abilities to process 
information. Here, entanglement is used as a resource which needs to be shared 
between several parties, eventually placed at remote locations. However en- 
tanglement is not the only manifestation of quantum correlations. Notably, also 
separable quantum states can be used as a shared resource for quantum commu- 
nication. The experiment presented in this paper highlights the quantumness 
of correlations in separable mixed states and the role of classical information 
in quantum communication by demonstrating entanglement distribution using 
merely a separable ancilla mode. 

The role of entanglement in quantum information is nowadays vividly demonstrated in 
a number of experiments. A pair of entangled quantum systems shared by two observers 
enables to teleport [8] quantum states between them with a fidelity beyond the boundary 
set by classical physics. Concatenated teleportations ^j can further span entanglement over 
large distances [10J which can be subsequently used for secure communication [TT]. An a 
priori shared entanglement also allows to double the rate at which information can be sent 
through a quantum channel [12] or one can fuse bipartite entanglement into larger entangled 
cluster states being a 'hardware' for quantum computing [13]. 

The common feature of all entangling methods used so far is that entanglement is either 
produced by some global operation on the systems that are to be entangled or it results from 
a direct transmission of entanglement (possibly mediated by a third system) between the 
systems. Even entanglement swapping [HIE], capable of establishing entanglement between 
the systems that do not have a common past, is not an exception to the rule because also 
here entanglement is directly transmitted between the participants. 

However, quantum mechanics admits conceptually different means of establishing entan- 
glement which are free of transmission of entanglement. Remarkably, creation of entangle- 



ment between two observers can be disassembled into local operations and the communica- 
tion of a separable quantum system between them [JS]. The impossibility of entanglement 
creation by LOCC is not violated because communication of a quantum system is involved. 
The corresponding protocol exists only in a mixed-state scenario and obviously utilizes less 
quantum resources in comparison with the previous cases because communication of only a 
discordant [T6HT8] separable quantum system is required. 

In this paper, we experimentally demonstrate the entanglement distribution by a sepa- 
rable ancilla [T3] with Gaussian states of light modes [TH]. The protocol aims at entangling 
mode A' which is in possession of a sender Alice with mode B' held by a distant receiver 
Bob by local operations and transmission of a separable mediating mode C from Alice to 
Bob (see Fig. [I]). This requires the parties to initiate their modes to a specific correlated 
but fully separable Gaussian state. Once the resource state Pabc is established no further 
classical communication is needed to accomplish the protocol. To emphasize this, we at- 
tribute the state preparation process to a separate party, David. Note, that this resource 
state preparation is performed by LOCC only. No global quantum operation with respect 
to David's separated boxes is executed at the initial stage and no entanglement is present. 

David's state Pabc has to be equipped with the property that by a local operation on 
subsystem (AC) it can be transformed to a state Pa'bc which is separable with respect to 
mode C, in order to guarantee that no entanglement is transmitted from Alice to Bob. It 
is also separable with respect to mode B', because no entanglement is a priori shared by 
Alice and Bob and the operation cannot create it. Mode A', however, has to be entangled 
with the composite subsystem (BC) so that the transmission of mode C could establish 
entanglement between Alice and Bob. It means that despite the impossibility to entangle 
Alice with Bob by LOCC, we are requested to create entanglement of mode A' with a 
delocalized two- mode system (BC), where mode C belongs to Alice and mode B belongs to 
Bob by a local operation on Alice's part of the initial fully separable state. What is more, 
the entanglement of mode A' with two-mode system (BC) should be localizable between 
modes A 1 and B' by a local operation on Bob's side. The entanglement carried by the global 
state is in fact bound and the localization renders its activation [20J which de facto crowns 
the entanglement distribution. 



The protocol 

The task can indeed be accomplished in the three-step protocol [HI] depicted in Fig. [IJ 
We assign the preparation of the fully separable three-mode Gaussian resource state, 
with properties described above, to David. Initially modes A and C are prepared in a 
momentum squeezed and position squeezed vacuum state, respectively, with quadratures 
xa,c = e±r %Ac-> P A -P = eTr PACi whereas mode B is in a vacuum state with quadratures 
x b = %b an d Pb — Pb ■ Here r is the squeezing parameter and the superscript "(0)" de- 
notes the vacuum quadratures. David then exposes all the modes to suitably tailored local 
correlated displacements [2T] (see also caption of Fig. [l]) : 

Pa -> Pa- P, x c ->■ x c + x, 

x B -> x B + V2x, p B -)> Pb + V*2p, (1) 

where the uncorrelated classical displacements x and p obey a zero mean Gaussian distri- 
bution with the same variance (e 2r — l)/2. The state has been prepared by LOCC across 
A|S|C splitting and hence is fully separable. 

In the second step, modes A, C of the resource state are passed to Alice and mode B to 
Bob. Alice superimposes modes A and C on a balanced beam splitter BSac- The beam 
splitter cannot create entanglement with mode B and hence the state is separable with 
respect to B\A'C splitting. Moreover, the state also fulfills the positive partial transpose 
(PPT) criterion [22} [23] with respect to mode C and hence it is also separable across C'\A'B 
splitting [24j as required (see appendix). 

In the final step, Alice sends mode C to Bob who superimposes it with his mode B on 
another balanced beam splitter BSbc- The presence of the entanglement between modes 
A' and B' can be certified by the sufficient criterion for entanglement by Duan et al [25J, 
T = ((xa 1 — x b') 2 ) + {{Pa' + Pb 1 ) 2 ) < 2. In the present case the total variance T reads as 

T-i + (^l)V-D + (^l)V 

and it satisfies T < 2 if r < ln[(4 + 3y / 2)/2] = 1.42. It clearly demonstrates the successful 
distribution of entanglement for a broad range of squeezing parameters r. This also confirms 
that mode A' was entangled with the two-mode system (BC) before the second beam 




FIG. 1: Sketch of the entanglement distribution protocol with continuous variables. 

David prepares a momentum squeezed vacuum mode A, a position squeezed vacuum mode C and 
a vacuum mode B. He then applies random displacements (green boxes) of the x quadrature (hor- 
izontal arrow) and the p quadrature (vertical arrow) as in (fiT), which are correlated via classical 
communication channel (green line). David passes the modes A and C to Alice, and mode B to 
Bob. Alice superimposes the modes A and C on a balanced beam splitter BSac an d communicates 
the separable mode C to Bob. This is vital for the global entangling operation implemented by 
the protocol. Bob superimposes the received mode C with his mode B on another balanced beam 
splitter BSbc'i which establishes entanglement between the output modes A' and B' (black lemnis- 
cata). Note the position of the displacement on mode B. In the original protocol the displacement 
is performed before BSbc, which is depicted by the corresponding box with dashed green line. In 
this case, David provides Alice and Bob with a fully separable three-mode Gaussian resource state 
exhibiting the specific correlation properties described in the text. Sending of the separable ancilla 
mediates the generation of entanglement between the modes A' and B' . Equivalently, this displace- 
ment on mode B can be performed after BSbc (dashed arrow indicates the respective relocation 
of displacement) on mode B' , and even a posteriori on a measured mode B' , as implemented in the 
experiment. The recovered entanglement is shown by the black lemniscata and confirmed by the 
experimental results. The solid red and green lines highlight the dual quantum-classical commu- 
nication channel effectively present in the protocol: separable quantum communication (red) and 
the flow of classical information (green) , respectively. 



splitter, because otherwise the establishment of entanglement across A'\B' splitting would 
not be possible. 

In the experimental implementation, instead of quadratures we will work with polariza- 
tion variables described by Stokes observables (see e.g. [261 127]). We choose the state of 
polarization such that mean values of Si and S2 equal zero while (5*3) ^> 0. This configura- 
tion allows to identify the "dark" Si-jS^-plane with the quadrature phase space. Sg,Sg +n /2 
in this plane correspond to Si, S2 renormalized with respect to S3 ~ S3 and can be associ- 
ated with the effective quadratures x, p. We also use the modified version of the protocol 
indicated in Fig. [T] by the dashed arrow showing the alternative position of displacement 
in mode B: The random displacement applied by David can be performed after the beam 
splitter interaction of B and C, even a posteriori on the measured mode B'. This is tech- 
nically more convenient and emphasizes that the classical information is sufficient for the 
entanglement recovery after the interaction of both modes A and B with the ancilla mode 
C and C, respectively. 

Experiment 

The experimental realization is divided in three steps, state preparation, measurement, 
and data processing, explained in detail in the methods section. David prepares two identi- 
cally, polarization squeezed modes and adds noise in the form of random displacements to 
the squeezed observables. The Stokes observable Sg is identified with p and Sg +7T / 2 with x. 
These two Gaussian modes in a mixed state are provided to Alice and then interfered on 
the balanced beam splitter BSac with a relative phase of | to ensure equal intensities of 
the output modes A' and C. 

The modulation patterns applied to modes A and C are chosen such that the two-mode 
state pa'c is separable. The states involved are Gaussian quantum states and, hence, are 
completely characterized by their first moments and the covariance matrix 7 comprising all 
second moments (see appenix). The covariance matrix ja'C f° r this two-mode state has 



lA'C 



(2) 



been measured to be: 

/ 20.90 ±0.0087 1.102 ±0.0091 -7.796 ± 0.0069 -1.679 ±0.0076\ 

1.102 ±0.0091 25.30 ±0.013 1.000 ± 0.0071 14.63 ± 0.0091 

-7.796 ± 0.0069 1.000 ±0.0071 20.68 ± 0.0093 0.8010 ±0.011 

\-1.679± 0.0076 14.63 ±0.0091 0.8010 ± 0.011 24.65 ± 0.0073 / 

A necessary and sufficient condition for separability of a Gaussian state pxy of two modes 
X and Y with the covariance matrix 'Jxy is given by the positive partial transpose (PPT) 
criterion |22, 



,(T Y ) 



7>T+^2>0, 



a-, 




(3) 
(4) 



where Jxy * s the matrix corresponding to the partial transpose of the state pxy with respect 
to the mode Y (see appendix). Effects that could possibly lead to some non-Gaussianity of 
the utilized states are discussed in detail also in the appendix. The state described by Ja'C 
fulfils the condition (3j as the eigenvalues (39.84, 28.47, 13.85, 9.371) of \t^%! + ifi 2 ) are 
positive, hence mode C remains separable after BSac- Bob mixes the ancilla mode C with 
a vacuum mode B on another balanced beam splitter, and performs a measurement on the 
transmitted mode B' . The measured two-mode covariance matrix of the output state ja'B' 
is given by: 

/ 19.95 ±0.011 1.025 ±0.016 -4.758 ± 0.0050 -1.063 ± 0.005l\ 

1.025 ±0.016 22.92 ±0.012 0.9699 ± 0.0047 9.153 ± 0.0058 

-4.758 ± 0.0050 0.9699 ± 0.0047 9.925 ± 0.0048 0.2881 ± 0.0047 

\-1.063± 0.0051 9.153 ±0.0058 0.2881 ± 0.0047 11.65 ± 0.0038 / 

Again, the separability is proven by the PPT criterion (eigenvalues 28.24, 21.79, 8.646, 
5.756). Only as Bob receives the classical information about the modulation on the initial 
modes A and C from David, he is able to recover the entanglement. After an appropriate 
post-processing of his measured raw data, Bob verifies that the Duan's criterion is vio- 
lated 1251 
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Anorm(SA',0° • g + Sb',0°) • A norm (5'A',90° " 9 ~ Sb',90°) < 1- 
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FIG. 2: Graph proving the recovered entanglement of the two mode state pA'B 1 - 
The experimental values for the Duan criterion are depicted in dependence of the gain factor g. 
Due to the attenuation of the mode B by 50%, a gain factor about 0.5 yields a value smaller than 
1, i.e. below the limit for entanglement (red line). The best entanglement value of 0.6922 ± 0.0002 
is shown for g opt = 0.4235 ± 0.0005. The inset zooms into the interesting section around the 
minimum. The depicted estimated errors are so small because of the large amount of data taken. 



That proves the emergence of entanglement (Fig: pi). The used gain factor g considers the 
slightly different detector response and the intentional loss of 50 % at Bob's beam splitter. 
The best entanglement value of 0.6922 ± 0.0002 is shown for g opt = 0.4235 ± 0.0005. This is 
the only step of the protocol, where entanglement emerges, thus demonstrating the remark- 
able possibility to entangle the remote parties Alice and Bob by sending solely a separable 
auxiliary mode C . 



Discussion 

Role of dissipation. The performance of the protocol can be explained using the struc- 
ture of the displacements Q. These displacements effectively render a coupling to the 
environment so that quantum correlations of a particular structure are shared with the 
environment and redistributed within the global system allowing to recover entanglement 
between the distant modes A' and B' in the final step. That is, after imposing the correlated 
displacements, the three modes are not in a product state any more. This can be interpreted 
as dissipation to a common reservoir, which can even lead to the creation of entanglement 
as shown in [251125]. 

Role of discord. All non-product separable Gaussian states have non-zero Gaussian dis- 
cord, that is, exhibit quantum correlations, which are beyond entanglement [30]. In our 
scheme, from the resource state prepared by David and through the protocol, the state of the 
modes held by Alice and Bob has non-zero discord which grows with increasing magnitude 
of displacements. It is in perfect agreement with the recent findings that for entanglement 
distribution by separable states to work, discord between modes A and B should be initially 
present. The entanglement gain is bound by the discord across the partition A'B\C after A 
and C have interacted [161 E] • Note that this latter requirement is directly related to the 
particular properties of our global state after the interaction of modes A and C on Alice's 
beam splitter described in the introductory part. That is, the possibility to distribute entan- 
glement by a separable ancilla can be attributed to the fact, that the state Pabc prepared 
by David is a separable quantum state possessing discordant correlations. Furthermore, the 
state pa'bc contains discord and entanglement across A'\BC splitting and is separable and 
discordant across C'\A'B splitting as required by the protocol. 

Role of classical information. Entanglement distribution without sending entanglement 
also highlights vividly the important role played by classical information in quantum in- 
formation protocols. Classical information lies in our knowledge about all the correlated 
displacement involved. This allows the communicating parties (or David on their behalf) to 
adjust the displacements locally to recover through clever noise addition quantum resources 
initially present in the input quantum squeezed states. Mode C transmitted from Alice to 
Bob carries on top of the sub-shot noise quadrature of the input squeezed state the displace- 
ment noise which is anticorrelated with the displacement noise of Bob's mode. Therefore, 



when the modes are interfered on Bob's beam splitter, this noise partially cancels out in the 
output mode B' when the light quadratures of both modes add. Moreover, the residual noise 
in Bob's position (momentum) quadrature is correlated (anticorrelated) with the displace- 
ment noise in Alice's position (momentum) quadrature in mode A' , again initially squeezed. 
Due to this the total variance of two modes A' and B' drops below the value for separable 
states and thus entanglement between Alice's and Bob's modes emerges. The difference be- 
tween the two protocol versions lies merely in the way how classical information is used. In 
the original protocol, the classical information is retained by David and he is responsible for 
clever tailoring of correlated noise. Bob evokes the required noise cancellation by carrying 
out the final part of global operation via superimposing his mode with the ancilla on BSbc- 
In the experimentally implemented protocol, David shares part of his information with Bob, 
giving Bob a possibility to activate entanglement a posteriori, by using his part of classical 
information after the quantum operation is carried out (see Fig. [I]). Thus entanglement 
distribution in our case is truly performed via a dual classical and quantum channel, via 
classical information exchange in combination with the transmission of separable quantum 
states. 

Methods 

State Preparation The two identically squeezed modes A and C required by the protocol 
are implemented as two polarization squeezed modes [26j [3TT433] . Each of these modes 
is generated by launching two orthogonally polarized femtosecond pulses (~200fs) with 
balanced powers onto the two birefringent axes of a polarization maintaining fiber (FS- 
PM-7811, Thorlabs, 13m). The pump source is a soliton-laser emitting light at a center 
wavelength of 1559 nm and a repetition rate of 80 MHz. By exploiting the optical Kerr 
effect of the fibers, the orthogonally polarized pulses are individually quadrature squeezed 
and subsequently temporally overlapped with a relative phase of 7r/2, resulting in a circular 
polarized light beam. The relative phase is actively controlled using an interferometric 
birefringence compensator including a piezoelectric transducer (PZT) and a locking loop 
based on a 0.1 % tap-off signal after the fiber. 

The modulation of each polarization squeezed mode in the dark-plane is implemented 
using electro-optical modulators (EOM). By applying a sinusoidal voltage V mo d, the bire- 
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fringence of the EOMs changes at a frequency of 18.2 MHz. We adjust the half-wave plates 
in front of the EOMs such that the modulation occurs along the direction of the squeezed 
observable. An additional quarter-waveplate is necessary to compensate for the residual 
birefringence of the EOM, even if no modulation voltage is applied. 

Finally, the two identically prepared modes A and C are interfered on a balanced beam 
splitter (BSac) with a fixed relative phase of tt/2 by controlling the optical path length of 
one mode with a PZT and a locking loop. This results in equal intensities of both output 
modes. 

Measurement Process To study the correlations between mode A' and B', multiple 
pairs of Stokes observables (Sa',6, Sb\o) are measured. In particular, the covariance matrix 
jA'B' is determined by measuring five pairs of observables: (Sa',o°, Sb',o°), (Sa',9o°, Sb',o°), 
(Sa',o°, Sb>, 90°), (Sa',9o°, Sb', 90°) and (Sa>a5°, <Sb',45°)- Here, 9 is the angle in the Si-iS^-plane 
between Sq° and Sg. The same procedure can be applied for the determination of ^a'C- 
The uncertainties indicated for the measured matrices have been determined by dividing 
the data set into 10 distinct parts and calculating the covariance matrices for the individual 
parts. The standard deviation of the elements of the covariance matrices gives the estimated 
statistical error. 

For the measurements of the different Stokes observables, we use two Stokes measurement 
setups, each comprising a rotatable half-wave plate, a Wollaston prism and two balanced 
detectors. The difference signal of one pair of detectors gives one Stokes observable Sg in 
the (S'l-^-plane, depending on the orientation of the half-wave plate. 

The signals are electrically down-mixed using an electric local oscillator at 18.2 MHz, 
which is in phase with the modulation used in the state preparation step. With this detection 
scheme, the modulation translates to a displacement of the states in Si-^-plane. The 
difference signal is low pass filtered (1.9 MHz), amplified and then digitized using an analog- 
to-digital converter card (GaGe Compuscope 1610) at a sampling rate of 10M samples/s. 
After the measurement process we digitally low pass filter the data by an average filter with 
a window of 10 samples. 

Due to the ergodicity of the problem, we are able to create a Gaussian mixed state 
computationally from the data acquired as described above. By applying 80 different 
modulation depths V mo d to each of the EOMs we acquire a set of 6400 different modes. 
From these set of modes we take various amounts of samples, weighted by a two dimensional 
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Gaussian distribution (see also Fig. [5]). The post-processing for the recovery of the entan- 
glement is performed on the measured raw data of mode B' . Therefore, the displacement 
of the individual modes caused by the two modulators is calibrated. By means of this 
calibration, suitable displacements are applied digitally. The classical noise inherent in 
the mode B' is completely removed. A part of the classical noise associated with Sa',o° is 
subtracted from Sb',o°, while the same fraction of the noise in Sa>, 90° is added to Sb',9o°. In 
this way, the noise partially cancels out in the calculation of Duan's criterion and allows 
to reveal the entanglement. We chose the fraction as in (fTl), which is compatible with the 
separability of the transmitted mode C from the subsystem (A'B) in the scenario with 
modulation on mode B before the beam splitter BSbc- 
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Note 

Recently, an experiment has been presented in [33] , which is based on a similar protocol. 
The main difference consists in the fact that it starts with entanglement which is hidden 
and recovered with thermal states. For this implementation no knowledge about classical 
information has to be communicated to Bob, besides the used thermal state. By contrast 
the setup presented in this work exhibits entanglement only at the last step of the protocol. 
Thus both works give good insights on different aspects of the theoretically proposed 
protocol [T9J. Another independent demonstration of a similar protocol based on discrete 
variables was recently presented in [35] . 
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FIG. 3: Sketch of the experimental setup Used abbreviations: HWP: Half-wave plate; QWP: 
Quarter-wave plate; EOM: Electro-optical modulator; BS1/BS2: Beam splitter; WS: Wollaston 
prism; (State Preparation) The polarization of two polarization squeezed states (A and C) is 
modulated using EOMs and sinusoidal voltages (18.2 MHz) from a function generator (dotted 
lines). The HWPs before the EOMs are used to adjust the direction of modulation to the squeezed 
Stokes variable, whereas the QWPs compensate for the stationary birefringence of the EOMs. The 
such prepared modes interfere with a relative phase of tt/2 on a balanced beam splitter (BS1). In 
the last step of the protocol, the mode C" interferes with the vacuum mode Bona second balanced 
beam splitter (BS2). (Measurement Process) A rotatable HWP, followed by a WS and a pair of 
detectors from which the difference signal is taken allows to measure all possible Stokes observables 
in the S1-S2 -plane. To determine the two-mode covariance matrix ^a'B 1 an necessary combinations 
of Stokes observables are measured. Removing the second beam splitter of the state preparation 
allows us to measure the covariance matrix of the two- mode state (A'C). (Data Acquisition) To 
achieve displacements of the modes in the Si-S^-plane we electronically mix the Stokes signals with 
a phase matched electrical local oscillator (18.2 MHz). The resulting signals are low pass filtered 
(1.9 MHz) and amplified before they are sampled by an analog-to-digital converter with a sampling 
rate of 10 M samples/s. The first step of the digital post processing made on the computer is an 
average filter with a window of 10 samples to further low pass filter the signals. After that the 
Gaussian mixed state is prepared. 
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FIG. 4: Histogram plots for Sa',o° of the Gaussian mixed state (blue) and three exem- 
plary individual modes. The transparent blue histogram depicts the Gaussian distribution of 
the mixed state, it is normalized to its maximum value. The red curve is a Gaussian fit to the 
histogram. The other three plots show the distributions for three out of 6400 displaced individual 
modes, the mixed state is composed of. They visualize the fact that the displacement of the indi- 
vidual modes does not affect the shape of their distribution. Gaussian functions have been fitted 
to them (red curves). The normalization is chosen such that they can be depicted in the same plot 
as the histogram of the mixed state. 
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Appendix A: Gaussian states 

We implement the entanglement distribution protocol using optical modes which are sys- 
tems in infinitely-dimensional Hilbert state space. An A-mode system can be conveniently 
characterized by the quadrature operators Xj,pk, j,k — 1,2, ... ,N satisfying the canonical 
commutation rules [xj,pk] = i8jk which can be expressed in the compact form as 

[| i ,|fc]=ifi i j fc . (Al) 

Here we have introduced the vector of quadratures £ = (xi,pi, . . . , x^,Pn) and 

N ( l\ 

V N = J, J=l , (A2) 

is the symplectic matrix. 

The present protocol relies on Gaussian quantum states. As any standard Gaussian 
distribution, a Gaussian state p is fully characterized by the vector of its first moments 



d = Tr (# j , (A3) 

and by the covariance matrix 7 with elements 

ljk = Tr[p{^- - djt, i k - 41}] , (A4) 

where {A, B} = AB + BA is the anticommutator. A real symmetric positive definite matrix 
describes a covariance matrix of a physical quantum state if and only if it satisfies the 
condition [36] : 

7 + m > 0. (A5) 

The separability of Gaussian states can be tested using the positive partial transpose 
(PPT) criterion. A single mode j is separable from the remaining N — 1 modes if and only 
if the Gaussian state p has a positive partial transposition p Tj with respect to the mode j 
[2H E7j. On the level of the covariance matrices, the partial transposition is represented 
by a matrix Aj = (©j-^ =1 ljj © a z , where a z = diag(l,— 1) is the diagonal Pauli z- 
matrix of mode j. The matrix 7™ corresponding to a partially transposed state p Tj reads 
ry(Tj) — Aj'jAj. In terms of the covariance matrix, one can then express the PPT criterion in 
the following form. A mode j is separable from the remaining modes if and only if [2U EZ] 

j(T d ) + [Q > 0. (A6) 
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The PPT criterion ( A6 ) is a sufficient condition for separability only under the assumption 
of Gaussianity of all the states involved. In our experiment, however, non-Gaussian states 
can be generated for which this criterion represents only a necessary condition. Therefore it 
can fail in detecting entanglement. 

Appendix B: Analysis of non-Gaussianity 

There are two sources of imperfections in our experimental set up that are potential 
sources of non-Gaussianity. These are phase fluctuations and the modulation of the initial 
squeezed states before the first beam splitter. They are discussed in the following sections. 

1. Phase fluctuations 

The experiment includes an interference of the modes A and C on a beam splitter, which 
is the first beam splitter BSac m the protocol. Imperfect phase locking at this beam splitter 
might cause a phase drift resulting in a non-Gaussian character of the state pA'c after the 
beam splitter. The phase fluctuations can be modelled by a random phase shift of mode 
A before the beam splitter described by a Gaussian distribution P{4>) with zero mean and 
variance a 2 . Denoting the operator corresponding to a beam splitter transformation as U 
and the phase shift on mode A as Va{4>), the state pA'c can be linked to the state Pac 
before the onset of phase fluctuations as 

Pa>c> = / P(0)W^ 4 (0)p AC V^(0)W t rf0. (Bl) 



Hence we can express the measured covariance matrix ^a'c given in Eq.(2) of the main letter, 
and the vector of the first moments d! of the state pA'c i n terms of the covariance matrix 
7 ac and the vector of the first moments d of the input state pac- For this it is convenient 
to define matrices D and D' of the first moments with elements Dij = didj and D[, = d^d'j, 



i, j — 1, . . . ,4. Using Eq. (Bl) and after some algebra, one gets the transformation rule for 
the matrix of the first moments in the form D' = £/X-D£[/ T , where U describes the beam 

CT 2 CT 2 

splitter on the level of covariance matrices. S = diag(e~~, e~~ , 1, 1) is a diagonal matrix. 



Similarly we get the covariance matrix 

lA!C > = U (S 7 acS + 7T © 0) U T , (B2) 
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where is the 2x2 zero matrix and 

(^ _ p -° 2 y 1 - p~ 2a2 

n = ± '-{A + a) + J(A + a)J T . (B3) 

Here the matrix A is the 2x2 matrix with elements Ay = ( r YAc)ij, hj = 1,2, a is the 2x2 



matrix with elements ocm = 2Dij, i,j = 1, 2, and J is defined in Eq.(A2). Similar to Ref, 



we can now invert the relation (B2) and express the input covariance matrix jac y i a the 



output covariance matrix ja'C an d the first moments after the beam splitter BSac as 



lAC = Y>- l U T lA ,c>UY>- 1 + 7f © 0, (B4) 



where 



7f = ^ ^~i A + «) + ^ J (^ + «) jT - ( B5 ) 

The 2x2 matrices A and a possess the elements Aij = {U T ^A'C'U)ij and 6>ij = 2{U T D'U)ij, 
i,3 = l,2. 

Our estimate for the variance of the phase fluctuations is a 2 = 0.02° and the vector d! of 
the measured mean values of the state Pa'c reads 

d! = (-0.208, 9.876, 13.32, 1.78). (B6) 



By substituting these experimental values for a and d' in Eq. (B4) and using the beam 



splitter with the measured transmissivity T = 0.49 we get a legitimate covariance matrix 



7 Ac before the phase fluctuations as can be easily verified by checking the condition (A5). 
Provided that the state with the covariance matrix 'Jac is classical it can be expressed 
as a convex mixture of products of coherent states. Gaussian distributed phase fluctuations 
and a beam splitter preserve the structure of the state, hence the state after the first beam 
splitter cannot be entangled. The covariance matrix / jac determines a physical Gaussian 
quantum state. Moreover, the covariance matrix possesses all eigenvalues greater than one 
and therefore the state is not squeezed [SB] which is in a full agreement with the fact that 
modulations of modes A and C completely destroy the squeezing. It then follows that this 
Gaussian state is classical and it therefore transforms to a separable state after the first 
beam splitter. 
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Inversion (B4) thus allows us to associate a Gaussian state before the phase fluctua- 



tions with the covariance matrix 'Ja'c measured after the first beam splitter. The sepa- 
rability properties of the state after the beam splitter can then be determined from the 
non-classicality properties of this Gaussian state. 

2. Gaussianity of the utilized states 



We have paid great attention on the modulations on modes A and C to preserve Gaussian 
character of the state pA'C- Our success can be visually inspected at the examples in Fig. 4 
of the main letter, which illustrates that both the modulation and the subsequent Gaussian 
mixing faithfully samples the required Gaussian shape. 

Besides this raw visual check we have also tested quantitatively Gaussianity of the in- 
volved states by measuring higher-order moments of the Stokes measurements on modes A' 
and C . Specifically, we have focused on the determination of the shape measures called 
skewness S and kurtosis K defined for a random variable x as the following third and fourth 
standardized moments 

(B7) 



S=^, K=^. 



where /^ = ((x — (x)) ) is the fcth central moment, (x) is the mean value and s = \fji2 is 
the standard deviation. 

Skewness characterizes the orientation and the amount of skew of a given distribution and 
therefore informs us about its asymmetry in the horizontal direction. Gaussian distributions 
possess skewness of zero. The exemplary values of skewness for various measurement settings 
are summarized in the Table IB 

TABLE I: Skewness S for Stokes measurements on modes A' and C in different measurement 
directions. 



Measurement 



Sa' 



o° 



SkewnessxlO 3 6.240 ±0.781 



Sa>, 



90° 



S C ' 



o° 



Sc, 



90° 



-1.478 ± 0.563 10.123 ± 0.727 1.106 ± 0.830 



The skewness can vanish also for the other symmetrical distributions, which may, how- 
ever, differ from a Gaussian distribution in the peak profile and the weight of tails. These 
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differences can be captured by the kurtosis which is equal to 3 for Gaussian distributions. 
The exemplary values of kurtosis for various measurement settings are summarized in the 
Table [TT1 

TABLE II: Kurtosis K for Stokes measurements on modes A' and C in different measurement 
directions. 



Measurement 


Sa',o° 


Sa',90° 


Sc,o o 


Sc, 90° 


Kurtosis 


2.971 ±2.211 x 1(T 3 


2.986 ± 1.852 x 1(T 3 


2.972 ± 1.978 x 1CT 3 


2.992 ± 1.568 x 10~ 3 



The tables reveal that the measured probability distributions satisfy within the experimental 
error the necessary Gaussianity conditions S — and K = 3. More sophisticated normality 
tests can be performed, which is beyond the scope of the present manuscript. 
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